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Abstract. We present two years of monitoring observations of the extremely variable quasar J1819+3845. We 
observe large yearly changes in the timescale of the variations (from ~ 1 hour to ~ 10 hours at 5 GHz). This annual 
effect can only be explained if the variations are caused by a propagation effect, and thus affected by the Earth's 
relative speed through the projected intensity pattern. To account for this effect, the scattering plasma must have 
a transverse velocity with respect to the local standard of rest. The velocity calculated from these observations is 
in good agreement with that obtained from a two telescope delay experiment (Dennett-Thorpe & de Bruyn 2001). 
We also show that either the source itself is elongated, or that the scattering plasma is anisotropic, with an axial 
ratio of >6:1. As the source is extended on scales relevant to the scattering phenomenon, it seems plausible that 
the anisotropy is due to the source itself, but this remains to be investigated. From the scintillation characteristics 
we find that the scattering material is a very strong, thin scatterer within ~ten parsecs. We determine a source 
size at 5 GHz of 100 to 900 microarcsecs, and associated brightness temperatures of 10^° to lO^^K. 
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1. Introduction 



' Variability in radio sources can be due to either intrin- 
, sic variations in the source itself, or due to a propagation 
effect known as scintillation. Disentangling these two ef- 
,, fects is critical to our understanding of the central regions 
of quasars and radio galaxies, as in both cases the vari- 
ability gives dimensions of the emitting region, but the 
interpretation (and the implied size) are quite different. 

Intrinsically variable extragalactic radio sources in- 
clude blazars, which exhibit changes in flux density of 
up to an order of magnit ude or more occurring over 
decades l|Aller et al.l Il985l). and smaller outbursts on 
short er timescales (e.g. IPohl et all Il996t IStevens et all 
Il995|l . Intraday variability is a phenomenon at GHz fre- 
quencies with t ypically a few percen t variations on day 
times cales (e.g. iHeeschen et all Il987t lOuirrenbach et all 
|2000|) . An intrinsic explanation requires the sources to 
have brightness temperatures far in excess of lO^^K ~ up 
to 10" - lO^^K. 

Scintillation occurs due to changes in electron density 
in the intervening material. Interplanetary scintillation is 
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caused by the solar wind. The variations of the flux density 
of quasars close to the Sun was used to deduce that some 
quasars had components of arcsecond size jHc wish et alJ 
Il964l ICohen et alJll967t iReadhead fc HewishI Il97fil) - a 
claim since confirmed by radio interferometry. Interstellar 
scintillation is caused by the material throughout the 
Galaxy, and is respon sible both for many effects seen in 
pulsars l)Rickettl[l990l and references therein) and low 
frequenc y variability (L FV) of some extragalactic radio 
sources l|Hunsteadll97^ . The compact structure in LFVs 
is expected to be of the order of a fe w milliarcsec, and is in - 
deed seen in VLBI observations(e.g. rSpangler et al.lll99^ . 
The scintillation regimes generally appropriate for extra- 
galactic sources are the (broad-band) weak and refractive 
regimes. (Pulsar observations are generally concerned with 
narrow-band diffractive phenomena. See Narayan, 1992 
for a review of scintillation). 
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Although refractive ISS is now commonly believed to 
be a major contributor to the variability at low frequen- 
cies dSoangler et al.ll989lll99ilBondi et alJl996l) . as first 
suggested by iRickett et al.l l)l984) . its relevance to ob- 
served variation s at GHz frequencies ha s been the subject 
of some debate llWagner fc Witzellll99,^ . in particular be- 
cause it cannot account for the reported correlation be- 
tween optical and radio variability, which have been ob- 
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serve d in a few sources l| Wagner et all Il996t IPeng et all 
|2000|) . The origin of the IDV in B 0917-^624, one of the 
most variable classical IDV sources, with 10-15% varia- 
tions, has been discussed by a variety of authors with dif- 
ferent views (e.glOuirrenbach et_al.ll_989HQian et ah 1991 



Simonettilll9 9l': Sha pirovskavalll99lt iRickett et al.lll995 : 
Kraus et al Tll99a) . An intrinsic explanation for the vari- 



ability in this source gives an apparent brightness temper- 
ature of lO^^K, which can be brought under the Compton 
limit by a Doppler factor of ~ 17 ((Ouirrenbach et all 
Il989l) . 

PKS 0405-385 was the first extragalactic radio 
source to show variations at GHz frequencies that 
were fairly unambiguously int erpreted as due to ISS 
I Kedziora-Chudczer et alJ 19971) . Using the simplest as- 
sumptions about the location and the velocity of the scat- 
terer, Kedziora-Chudzner et al. calculated the brightness 
temperature of PKS 0405-385 to be 5 x lO^^^K, or 1000 
times the inverse Compton limit. However, this figure de- 
pends critically on the assumptions about the scatterer. 

In this paper we investigate the properties of the flux 
density variations of the quasar J1819-I-3845. The varia- 
tions in this source are the most extreme known in the 
radio sky. We previously interpreted the variability as due 
to interstellar scintillation (ISS), based on the extreme 
nature of the variations and their f requency dependence 
l|Dennett-Thorpe fc de Bruynll2000l hereafter Paper I). 
We have since proved them to be due to ISS by detect- 
ing a difference in the arrival times (of up to 100 sec) of 
the flux density variations at two wi dely spaced telescopes 
( Dennett-Thorne fc de Bruvnll200ll hereafter Paper II) . 

In Paper I we calculated the location of the scatterer 
using simple formulae and the observed timescale, depth 
of the intensity variations and the frequency of maxi- 
mum variations. We assumed that the scatterer was at 
rest w.r.t. the LSR, and thereby calculated in a straight- 
forward manner that the brightness temperature of the 
scintillating part of J1819-I-3845 was > 5 x lO^^K. This is 
an order of magnitude higher than the Compton limit- a 
fact which could be explained by Doppler boosting, but 
this expects rapid energy loss at higher radio frequencies. 
Furthermore all known sources with Doppler factors of 10 
or more are themselves intrinsically variable. 

We conducted a monitoring campaign on J1819+3845 
to determine the longevity and stability of the source and 
the velocity of the scatterer (an unknown in the T^ cal- 
culation) . Since the programme started we have measured 
the velocity of the scattering plasma by another means 
(Paper II) . With the monitoring campaign we also wished 
to try to 'resolve' the source by observing a change in 
scintillation characteristics throughout the year, as we 
sampled the source along different axes. Here we report 
on the monitoring campaign at 5 GHz and will use the 
conclusions in later papers which will discuss the (multi- 
frequency) structure of the source. 

In Sect. 121 we describe the observations (and include an 
optical spectrum). We give an overview of the results in 
Sect. 01 before undertaking an analysis of the light-curves 



(Sect.^. Sect.OandElgive the interpretation of the results 
by fitting models for plasma velocity and anisotropy to the 
data. In Sect. 7 and 8 give give some further consideration 
to the scattering medium and source stability. In Sect. 9 
we derive the parameters of the scattering medium, and 
the size of the source, taking into consideration effects of 
source resolution. We end with a summary of the results 
and concluding remarks. 

2. Observations 

2.1. Radio 

The observations were conducted with the Westerbork 
Synthesis Radio Telescope (WSRT) in the period May 

1999 to July 2001. The details of the observations are 
given in Table ^ The telescope generally observed in full 
array mode, with all antennae working at the same fre- 
quency, but on a number of occasions two frequencies were 
observed separately in split-array mode. From December 

2000 onwards most observations were taken in split-array 
configurations with 5 telescopes at 2.3 GHz, 4 telescopes at 
4.9 GHz and 5 telescopes at 8.5 GHz. The number of tele- 
scopes used at 4.9 GHz is indicated in Gol. 4 of Table ^ 
(Data at the other frequencies will be presented and anal- 
ysed as part of a forthcoming paper.) On-line integration 
times were set to 10 seconds. The total recorded band- 
width was 80 MHz, in 8 contiguous 10 MHz bands. In the 
split-array observations software constraints limited the 
number of independent bands to 6x10 MHz. 

Most observations consisted of long, uninterrupted 
tracks, of up to 12h duration. Each long track was pre- 
ceded and followed by short observations of polarized and 
unpolarized calibrators (generally 3C286/CTD93 before 
and 3C48/3C147 afterwards). Absolute flux density cal- 
ibration was performed assuming a flux density of 7.42 
Jy for 3G286 and 1.70 Jy for GTD93 at a frequency of 
4.874 GHz. The polarisation characteristics will be dealt 
with in a later paper. 

The reduction package newstar (Noordam et al., 
1994) was used to calibrate and examine the data. Bad 
bands, telescopes and baselines were excised where neces- 
sary. The amplitude and phase calibration solutions deter- 
mined for the calibrator were applied to J1819-f3845. The 
amplitude stability of the WSRT is generally very good. 
The errors generally contain three components which are 
due to pointing (< 1% at 4.9 GHz), system temperature 
and atmospheric opacity variations. The main observed 
variations in Tgys are due to variable spillover losses and 
atmospheric emission, both of which vary only slowly. 
The WSRT continuously measures Tgyg from stable noise 
sources, allowing these slow variations to be accurately re- 
moved, resulting in an overall array amplitude stability of 
better than 1% over periods in excess of 12 hours. 

The largest remaining uncertainty is due to the atmo- 
spheric opacity variations associated with (heavy) rainfall. 
The emission effects on the Tgys are taken care of but the 
associated opacity effects are not. The uncertainty in the 
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flux density scale due to opacity effects is less than 1%. 
Errors associated with heavy rainfall, occurring less than 
1% of the time, can, however, lead to an underestimate of 
the flux density of a few percent. 

Phase self-calibration was applied to the observations 
of J1819-f 3845 , with a 5 minute averaging interval. The 
source model of J1819-I-3845 at 4.9 GHz was a single point 
source at the phase centre, as there are no sources detected 
above 1 mJy within the primary beam. Previous observa- 
tions with the VLA had already confirmed that there is 
no structure in the source on scales in between 0.3 and 
3", the resolution of the WSRT at 4.9 GHz. The source 
flux density was then calculated from the real part of the 
complex visibility, averaged over the desired sampling in- 
terval (30 or 120 s). This requires that the phase is stable 
over the period of the self-calibration solution, a condi- 
tion seen to be held from inspection of the slowly chang- 
ing phase solutions. This method has the advantage that 
it does not require a noise bias correction which is nec- 
essary if scalar averaging would be used. Rare occasions 
of rapid atmospherically induced phase variations, were 
detected by comparison with the scalar averaged data: on 
27/08/1999, 08/09/2000 and 13/10/2000, when the phase 
fluctuations result in a minor (a couple of percent) time- 
dependent variation in the calculated flux density. 

Typically between 70 and 80 baselines contribute to 
the final datum. With an rms noise of 20mJy per base- 
line, for 120 s averaging times in one 10 MHz band, we then 
achieved a typical noise in the source flux density of about 
0.8 mJy for eight 10 MHz bands and 120 s integration. This 
0.8 mJy is the thermal noise error. For the 4 telescope (6 
baseline), 60 MHz observations from December 2000 on- 
wards the error is about 4 times larger, i.e. 3 mJy. The 
final flux density error is the quadrature sum of a thermal 
error and the slowly varying gain error discussed above 
which we estimated at typically 1%. With a median flux 
density of the source of 200 mJy we can see that the over- 
all error is never more than about 1.5%. The smooth light 
curves, in periods of slow scintillations, show that our es- 
timations of the errors are realistic. 

For the 10 minute observations conducted in 
September 1999, and in two periods in September 
2000 and November 2000, the errors in the absolute flux 
density scale outweigh the errors from the fit. 

2.2. Optical 

A spectrum was obtained at the WHT on 17 July 2000 
in three 1800 sec exposures. The ISIS spectrograph was 
used to obtain a spectrum from 3000A to 8OOOA, with a 
resolution of 0.86 and 2.9 A/pixel and central wavelengths 
of 4495A and 6708A on the blue and red arms respec- 
tively. Standard reduction procedures were followed using 
the IRAF longslit software package. The combined spec- 
trum, with the blue arm boxcar smoothed over 9 pixels, 
is shown in Fig^ The fitted line properties are shown in 
Table H 



Table 1. Observational parameters 



Date day MJD duration number notes 

number (hrs) tels 



month 


day numbers 


duration 


telescopes 


1999-09 


244,245,246,247,248, 


lOmin 


14 




252,254,256,261,262 






2000-09 


244,245,246,247,248,249 


3min 


14 


2000-10 


300,301,302,303,306 


3min 


14 




307,308,309,311,312,313 







We calculate the redshift to be 0.533±0.001. The spec- 
trum is fairly typical for quasars, both in the equivalent 
widths of the lines and in the continuum emission (see e.g. 
Miller et al. 1992; Baker 1997; Brotherton et al. 2001). 
We note the relatively small Balmer decrement (H/3/H7), 
possibly indicating little line of sight obscuration by dust. 
The [OH] luminosity (lO'^'^W) is comparable to moder- 
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Fig. 1. WHT spectrum of J1819-I-3845, with identified lines indicated 



ately powerful FRII radio sources (radio luminosities at 
151MHz of 1025-1027w/Hz/str)tWillott et al. 1999). We 
have observed the optical spectrum of this source twice. 
On both ocassions the optical magnitude was similar, and 
the lines strong and broad. From these observations it does 
not seem that J1819+3845 is a typical OVV; rather a typ- 
ical quasar. 

Table 2. Fitted line properties of optical spectrum 
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3. Overview of radio observations and comments 
on individual epochs 

In Fig. 12 we present the light-curves over more than two 
years of observations. In these curves the averaging time 
is 120 s. For the analysis we used 30 sec integrations, an 
example of which is plotted in Fig. |3 The interrupted 
light curves in September 1999 are due to rapid switching 
between two frequencies about 100-200 MHz apart. 

The rapid variations of May 1999 have become slower 
by June and July of that year, although they still are 
fast enough to show complete modulations in a few hours. 
By the end of August the modulations have become so 
slow that the source remains at an almost constant flux 
density over 6 hours. The fact that the source is still vary- 
ing, but with a very long timescale can be seen in Fig. [3 
which shows the variations over August and September, 



including the short (~ 10 mins) observations taken over 
the ten days following the quasi-stationary observations 
on August 27 (day 239). The source was observed close to 
its minimum flux density (60 mJy) on the next observa- 
tion (September 1, day 244), and on the following day at 
a flux density typical of the peak excursions (220 mJy). 

The modulations remain slow until December, when 
they quickly speed up again. By the end of 1999 there are 
rapid intensity variations with peak to peak timescales of 
around one and a half hours. The variations remain fast, 
up to and including the observations of 15 May 2000, one 
year after the first twelve hour observations. It can be seen 
that the character and timescale of the observations is ex- 
tremely similar to the observations of one year previously. 

The observations during May 2000 to July 2001 
broadly replicate the observations of the previous year: 
showing a slowing down of the variations between June 
and December, before speeding up again. Note also the 
very rapid speeding up of the scintillations in December 
of both years: a factor of 2 within an interval of just 3 
weeks. 

We have analyzed variations across the 80 MHz bands 
and found them highly correlated. The correlation coeffi- 
cients of > 0.99 are consistent with completely correlated 
signals with thermal noise. At 5 GHz the scintillation is 
clearly broad-band (Fig.Q. 

4. Analysis of light curves 

4.1. Modulation index 

The modulation index is the normalized rms intensity 
variations, -\/S(S' — S)'^/S. In the asymptotic theoretical 
limit, the modulation index of broad-band variations of a 
point source (a source subtending an angle much smaller 
than the scattering angle) reaches unity at a 'critical fre- 
quency'. Above and below this frequency the modulation 
index drops. 
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Fig. 2. Two years of monitoring observations of J1819+3845 at 5GHz. The first page covers 11 months, the second 
page 15 months. The annual change in the timescale of the modulations is clearly apparent. The horizontal axis shows 
the Hour Angle of the source at the WSRT, so that each panel covers to 12 hours. 
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We have calculated the modulation index for each 
of the observations and present the results in Table O 
However, during the period of slow modulations the 'mod- 
ulation index' calculated during one run has little mean- 



ing, as typically only one scintle is seen per observation. 
We therefore also present in Table 13 a modulation index 
calculated over the entire period July to October for both 
years. 
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Fig. 3. The light curve at 30 sec resolution from 05 
February 2000. 
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Fig. 4. The light curves from 17 July 1999 at 120 sec res- 
olution showing all eight 10 MHz bands. The lowest light 
curve corresponds to a central frequency of 4.839 GHz. 
The other bands, each 10 MHz higher in frequency than 
the previous, are displayed with an offset interval of 
25mJy. 



We calculate the error on the modulation index from 
Monte Carlo simulations. We used the observations which 
display over eight strong intensity maxima, and selected 
portions of these runs, to simulate the epochs where the 
variations are slower and we have few scintles. We calcu- 
late both an estimate of the scatter on calculated mod- 
ulation indices, and an estimate of the underestimation 
that is caused by observing fewer scintles (^Naravan et al.i 




220 240 260 280 300 320 

UT day number 

Fig. 5. Observations in August through October in 1999 
(top panel) and 2000 (bottom panel) . In the 6-hour obser- 
vation on 27 August 1999 (day 239), the source appears 
to maintain a near-steady flux density: but the variations 
have not stopped, merely slowed down. 



llQSgtlPeshpande &: Nitvanandall990|) . For an observation 
of N high maxima ('scintles'), the scatter on the calculated 
modulation index is given by ~ 0.3m/-\/iV, and the under- 
estimation of the modulation index is typically around half 
this. The error on the calculated mean flux density S is 
found to be ~ S/6N, or 2%, whichever is the larger. 

For the period of short observations we have calculated 
the mean flux densities and modulation indices. These are 
shown at the end of Table O As the observed flux density 
is not normally distributed, we again estimate the error 
on the calculated quantities by use of Monte Carlo simu- 
lations. This time we drew a number of flux densities ran- 
domly from an observation with many variations. We con- 
clude that the error on the calculated modulation index, 
TO, after N observations, can be expressed as ^ 2m/5^/N. 
The error on the calculated mean S can be expressed as 

From this analysis we conclude that the modula- 
tion indices at 5 GHz are not statistically distinguishable 
throughout the observing period (Figl^l). They remain 
constant at ^0.40. 

4.2. Timescales of variation 

The timescale of the variations can be defined in a number 
of ways. The timescale used in Paper I was defined by use 
of the structure function, as this is the quantity easily 
related to the theoretical studies of scintillation. However 
this measure is impossible to compute with less than one 
observed scintle. With very few observed modulations any 
measure of the timescale will obviously be dominated by 
statistical error. Therefore we have used independent ways 
of estimating the timescale in an attempt to minimize the 
effect of this on our conclusions. We describe these three 
ways of calculating the timescale below, and present the 
results in Col. 8 to 11 of Table El 



8 



J. Dennett-Thorpe and A.G. de Bruyn: Annual modulation in the scattering of J1819+3845 




400 



600 800 1000 

MJD - 51000 



1200 



Fig. 6. Modulation index as a function of observing date. 
The mean flux densities used in this plot correspond to 
the mean measured during the respective observation. The 
two dotted crosses mark the modulation index calculated 
using all the observations in the slow period between July 
and October. Our estimation of the errors, given by the 
vertical bars, is explained in the text. 



For our regularly sampled light curves f(i), i=l,2...n, 
the first order structure function is estimated as 



D){k) 



1 



w{t)w{i + k)[f{t + k)-f{t)f 



where w(i)=l if the measurement exists and is otherwise 
zero (Simonetti et al. 1985). The structure function anal- 
ysis is straightforward for those observations with a large 
number of modulations. In this case it can be seen that 
the structure function has saturated (see e.g. Fig. [Tj), and 
we can calculate the tiniescale as twice the time at the 1/e 
of this maximum value (which is D^(max)= 2m^). This is 
given by t in Col. 8 of Table El 

Following the error analysis of Simonetti et al. (1985) 
the observed value Dy(k) exceeds the true value by an 
additional quantity 2cr^^, where agj is the measurement 
noise variance (assuming white noise with correlation 
timescale smaller than the smallest time interval mea- 
sured) . We thus assess the timescale correcting each value 
of with the relevant noise variance, but the correction 
on the calculated timescale is insignificantly small (~1 sec 
for the most rapid variations). The error on the calculated 
structure function Dj{k) from this contribution is negli- 
gible. From the Monte Carlo simulations, we find that the 
(one sigma) error on this timescale due to statistics of 
finite observing period is ~ 0.25t/^/N. 

For observations with few (<5) scintles we estimate 
a timescale, from the structure function by assuming a 
mean modulation index for all epochs. By considering the 
modulation indices obtained in May 1999, and December 



1999 through May 2000, we take a modulation index of 
m=0.4 for all epochs and estimate the timescale as the 
time at which the structure function (or its extrapolation) 
reaches 2m^/e. This modification is given by t' in Col. 9 
of Table El 

As a final, and most direct, approach, we also esti- 
mated the timescale, t*, from the average time between 
two minima (or maxima) during the observation run. For 
runs with a small number of scintles, t* was calculated 
as the average of twice the time between a minimum and 
maximum. We estimate errors on this process by consid- 
ering the runs with >5 complete modulations. We cal- 
culated the cr on the timescale as measured from each 
single complete modulation independently. We found that 
the error on the measured timescale can be estimated as 
at* = 0.3t* /y/N, where N is the complete number of mod- 
ulations observed during the run. We calculated t* in or- 
der to assess the reliability of our error estimates on runs 
with few scintles, and their possible propagation into the 
final results. The quantitative results obtained using t* 
are quoted in Sect. El but we note here that the results 
so obtained were similar to those obtained by the other 
estimations. 

For the observations of 9 March 2000, it is not clear 
how many minima occur during the run, due to gaps in the 
observations. Considering the rapid change of flux density 
on this date, it is plausible that extra minima also oc- 
cur during the unobserved periods. The timescale quoted 
is based on this assumption, the error used in the fits is 
based on the difference between this and the timescale cal- 
culated if only 4 minima occur during the observations. 
We can obtain lower limits to t* from the duration of 
the observations during which the light curves showed lit- 
tle variation: 1999-08-27 (360 min), 1999-09-21 (300 min), 
1999-11-28 (400 min), 2000-08-27 (400 min). This informa- 
tion is not included in the fitting. 

5. Detection of peculiar velocity of scattering 
plasma 

Fig. El shows the changing timescale (t) as a function of 
epoch throughout the period of observations. The yearly 
changes in the timescale can also be clearly seen directly 
in the data (Fig.EJ- 

A change in the timescale of the scintillations with a 
year periodicity must be due to the changing speed at 
which the Earth cuts through the scintillation pattern 
projected onto the Solar System, as we move around the 
Sun. The velocity of the scintillation pattern, as viewed 
by a stationary observer on Earth, is given by the rela- 
tive velocity of the Earth and scattering plasma projected 
onto a screen perpendicular to the direction of the source, 
I'Eaith — I'piasma- We notc that such a model predicts no 
change in the modulation index. 

For most extragalactic sources, the vector WEarth itself 
will provide a substantial yearly modulation in the pat- 
tern speed. This is in contrast to pulsars, which produce 
scintillation patterns whose observed motion is dominated 
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Table 3. Derived quantities: Col. 5 is the number of scintles observed, N. Cols. 6 & 7 list the modulation index and 
associated error. Cols. 8, 9, 11 show the derived timescales: using the structure function and derived modulation index, 
using m==0.4, and an estimation based on the average time between extrema, respectively. 
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Combined data from slow season 
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by the pulsar's own velocity of typically several hundreds 
km/s. One exception is PSRB2016+28, which has a trans- 
verse proper motio n velocity comparable to the Earth's. 
iGuDt a et al.l l)l994|) found a variation in the scintillation 
timescales of pulsars, which they attributed to an annual 
variation caused by the Earth's motion. Unlike most ex- 
tragalactic sources, J1819-I-3845 is not expected to show 
a large annual modulation in the pattern speed. This is 



illustrated by the dotted line in Fig.|Hl J1819-f 3845 is just 
9° from the direction of Sun's motion with respect to the 
LSR, and 28° from the ecliptic pole, so the modulation 
due to the Earth's motion is slight. 

In order to fit the large change in timescale during 
the year we require that the scattering plasma itself has 
a velocity with respect to the solar rest frame (taken as 
moving at 19.7 km/s towards 18 07 50, -^30 00 52 (J2000)). 
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m = 0.42 
10.0 < t(char) < 12.0 
D1 slope ~ 1 .89 



time (mins) 

Fig. 7. A saturated structure function with observations 
from 18/03/2000. Crossed circles are points obtained di- 
rectly from the data, and the crosses are after correction 
for thermal noise (see text). The error on the calculated 
timescale is dominated by statistical error from a finite 
number of observed variations, not the thermal noise cor- 
rection. The solid lines indicate slopes cx t'^ and t (see 
Sect. 7). 



The change in timescale throughout the year is assumed 
to be entirely due to changes in our velocity through an 
isotropic scintillation pattern. Thus we assume that the 
size of the scintles, angular size of the source, and dis- 
tance to the scattering screen are constant throughout the 
years. However, we neither fix nor determine these values, 
as the change in timescale only depends on the changing 
projected velocity. 

We probed a region of parameter space for the pro- 
jected velocity of the scattering plasma between -50 km/s 
and +50 km/s in both RA and declination, and calcu- 
lated the best fit. The best fits were clearly constrained in 
this space, and the best weighted fit to the calculated 
timescales is given by the solid line in Fig. |H1 It shows that 
the change in timescale can be qualitatively explained by 
a velocity of the scattering plasma of 30kms~^. The 
model allows the timescale to be normalised according to 
the best- fit. The ellipses in Fig. [SJa correspond to the pro- 
jected velocities, for the cases of no peculiar plasma veloc- 
ity (dotted line) and a plasma velocity corresponding to 
the fit shown (solid line). The velocity vectors (using the 
fitted plasma velocity) are also shown for the beginning of 
each month. 

Although the model qualitatively describes the data, 
with an annual period of long modulations, it is clear that 
the fit is poor (reduced — 9-1 for goodness of fit 
parameter Q=10~'^^). The fit is similar using the other 
timescale estimates. The period of very slow modulations 




100 200 
\y day number 

Fig. 8. The timescales, t, (logarithmic scale) of the vari- 
ations as a function of observation date. The dotted line 
corresponds to the prediction if the scattering plasma has 
no peculiar velocity. The solid line gives the best fit to 
the data (using timescale estimate t) if the plasma is al- 
lowed a peculiar velocity. The lower panel (b) shows the 
data folded on a year. The annual modulation, as well 
as the width of the epoch of slow modulations can be 
clearly seen. The corresponding monthly velocity vectors 
t^Eaith — Wpiasma are shown in the ellipses: near-circular 
around the origin for no plasma motion, and offset for the 
fit shown (see also Fig [TT)l . 



extends for a much longer time than is allowed by this 
simple model. Nor is the sharp reduction in timescale in 
December, or sharp increase in timescale in July/ August 
each year replicated by this model. Furthermore the fitted 
plasma velocity does not agree well with that measured 
from a delay experiment (Paper II). However, a peculiar 
velocity of the scattering plasma is the only possible ex- 
planation for this annual variation. In the next section we 
consider further refinements and additions to the model 
which can explain the long period over which slower mod- 
ulations are observed. 



J. Dennett- Thorpe and A.G. de Bruyn: Annual modulation in the scattering of J1819+3845 



11 



6. Detection of source structure/ anisotropic 
scattering medium 

In the preceding section we assumed that the scintilla- 
tion pattern was isotropic. Anisotropic scintles could be 
produced by either an extended anisotropic source or an 
anisotropic scattering medium. An anisotropic scintilla- 
tion pattern would leave a distinctive trace in the form 
of an annual variation of the timescale. This would hap- 
pen as the Earth moves around the Sun, and we cut the 
pattern at varying angles through the anisotropic pattern, 
thereby introducing an annual variation in the timescale 
of modulations. 

In the absence of any peculiar velocity of the scatter- 
ing plasma, the Earth would cut through the scintillation 
pattern of J1819-f3845 at the same angle, but moving 
in different directions, approximately every six months. 
Thus, for this source in the absence of any peculiar scatter- 
ing plasma, the lengthening of the timescale would occur 
about every six months. It is therefore clear that source 
structure cannot be the only explanation of the observed 
changes in timescale. 

With the addition of a peculiar velocity of the scat- 
tering medium, however, the periods of slow modulations 
no longer occur every six months. It is therefore possible 
that both a peculiar plasma velocity and source struc- 
ture combine to produce the observed annual variation in 
timescale. 

We modelled an anisotropy in the scintillation pattern 
as an ellipsoid, introducing two new parameters: the po- 
sition angle and axial ratio of the scintles. Using these 
extra parameters we again performed a minimization 
over a wide range of parameter space. The resulting best 
fits are significantly improved over the fits without the 
anisotropy term (reduced is 1-5 for t; 5.2 for t*; 7.1 for 
t'). Fig.|51shows the probability contours. Although the fit 
was improved, good solutions were available over a range 
of parameters. It was found that the position angle and 
Vdcc were relatively well-constrained, but the axial ratio 
and vra were poorly constrained. 

In Paper II we directly measured the velocity of the 
scattering plasma using the time delay of signals arriv- 
ing at two telescopes: vra = —32.5 ± 0.5km/s; Vdoc = 
-1-15.5 ± 1 km/s.^ This result lies within the narrow re- 
gion of permissible VRA-Vdec space determined from these 
yearly monitoring observations (see Fig Ell. Constraining 
the velocity to be within the 3-cr errors of this measure- 
ment, we re-ran the fit and determined a probability pro- 
file for the shape and position of the ellipse. This is shown 
in Fig.llOb. and the results are listed in Table0] The best 
fit (at vra ~ —33.5 km/s; Vdcc = -t-13.5 km/s) has axial 
ratio 147^ N at position angle 83° ± 4° (N through E). 



^ Note added in proof: We have since realised that full con- 
sideration of anisotropy in the two-telescope experiment is also 
needed. The effect of this is that the velocity is not as well 
constrained as given here, but the main conclusions are not 
significantly altered. 



Use of the timescales obtained from the peak separation 
gives a similar result, although the formal fit is very poor. 

The introduction of two more parameters has clearly 
improved the fit sufficiently to justify their use, as can be 
directly seen in Fig. EH We also note that the dates on 
which the observations display very little variation over 
6 hours, and which therefore had no calculated timescale 
and were not included in the fit, occur during the periods 
of predicted slow modulations. 

The fit is illustrated in an illuminating way in Fig. IT21 
which shows the data compressed by the timescales pre- 
dicted from the model, such that the effects of plasma 
motion and anisotropy are removed. After this timescale 
adjustment, the similarity of the modulations throughout 
the two years can be clearly seen. The timescale of the vari- 
ations now looks similar throughout the period, indicating 
the model fit was suitable. The nature of the variations is, 
by eye, also similar throughout the observing period: there 
is no indication of big changes in the 'jaggedness', or the 
number and frequency of big excursions. What is clear 
from the figure however, is the increasing flux density of 
the brightest peaks in the light curves. The minimum flux 
densities seem to have changed little, showing at most a 
slight increase. This information is also reflected in Table 
El [Cols. 3 & 4]. As was shown in Fig. El there has been 
no change of rms variations over time. The change in ap- 
pearance of the light curves over time, in particular the 
increase in high maxima, is due to the increasing flux den- 
sity of the source (see Sect.jSJ. 

One manner of producing the ellipsoidal scintles is by 
an ellipsoidal source, as the scintillation intensity pattern 
produced by an extended source is the convolution of the 
scintillation pattern of a point source with the bright- 
ness distribution of the source itself. In the strong, re- 
fractive regime a resolved source (source size 9s > ^scatt) 
is observed to have a timescale of modulation larger than 
that of a point source as t oc 9s /Oscatt (Narayan, 1992), 
where Os is the width of the source. If we assume that the 
anisotropy is then due to an ellipsoidal source undergoing 
strong scattering, we may equate the ellipsoid calculated 
above with that of the source itself. 

We note that we are observing near the critical fre- 
quency where the asymptotic theory will not hold to great 
accuracy. However, as the asymptotic theory from the 
other side of the critical frequency, in the weak scattering 
regime, predicts r cx ^5^^, we take r cx ^5 as a reason- 
able approximation. The second complication is that to 
apply the theory as simply as we have done at least re- 
quires the source to be larger than the scattering disk in 
all directions. We believe however this is a reasonable as- 
sumption, as a source that was significantly smaller than 
the scattering disk in one direction should be expected 
to display different characteristics of scintillations as the 
velocity vector cuts across its minor axis. At this time it 
should display more of the second branch, or diffractive, 
scintillations. There is no clear evidence for such a change 
in character, such as a second, faster, timescale of modu- 
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Fig. 9. The probability distributions based on ^^-minimization fits, using t. For each plot the probabilities are cal- 
culated over the ranges shown in the other plot. The contours correspond to 68.3, 95.4, 99.73 and 99.99% certainties 
for a normal distribution of errors. The vertical line in the first panel represents the velocity and its errors derived in 
Paper II. Note that in these plots PA is defined as E through N 



lation. Thus for the anisotropy caused by an anisotropic 
source, in so far the choice of ellipsoidal model is a good 
representation, the values obtained (Tabled are expected 
to be close to the anisotropy of the source itself. 

An anisotropy in the scatterer could produce an ef- 
fect similar to the one described above for an anisotropic 
source. Anisotropy in the scattering medium is generally 
expected, and has been diagnosed by elongated scatter- 
broadened images, with axial ratios of up to 3 found 
through strong l y scattering regions of the Galaxy (e.g. 
Lo et al.l Il993t IPrail et all Il994t IWilkinson et"aLlll99^ 
^otteret al.l Il998[). and axial ratios > 5 in th e solar 
wind (e.e. iNaravan et"aLl Il989l: lArmstrong et al.lir99Ql . 
Compressible MHD turbulent models have been devel- 
oped in which the magnetic field gives rise to an elonga- 
tion of the electron-density fluctuations, with potentially 



very large axial ratios l|Lithwick fc GoldreichllioOlf) . The 
effect of this anisotropic medium on the scintillation ob- 
servables is greatest when the screen is thin, the mag- 
netic field perpendicular to the line of sight and with lit- 
tle change of magnetic field direction in the scatterer. In 
order to produce the > 6:1 axial rati o which would be re- 
quired for J1819-I-3845 , we see from lChandran fc Backer! 
I|2002|) that the magnetic field in the plane of the scattering 
screen must rotate through less than 7r/4 radians through 
the screen. This implies that if anisotropy in the plasma 
density structure causes the observed anisotropic scintles, 
the magnetic field must be well ordered through the scat- 
tering screen (i.e. along line of sight). The apparent >6:1 
axial ratio required is somewhat unusual for high galactic 
latitudes, but may be related to the unusual nature, and 
thinness of the screen. 
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Fig. 10. Fitting constrained by the plasma velocity (and error) from Paper II. Timescale estimate t has been used, (a) 
The probability distributions based on fits, showing the constraints placed on the axial ratio and position angle of 
the anisotropy of the scintles. (b) The resulting best fit to the timescales with epoch. 



Table 4. Results of fitting, using vra — 
Vdcc — 15.5 ± 3 km/s. Errors are 3-cr. 
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7. Grand Structure Functions 

In order to attempt to distinguish between the possibil- 
ities of source or plasma anisotropy, to get clues about 
the depth of the scattering plasma, we constructed grand 
structure functions {D^ (r)) . For this purpose we took both 
the mean calculated from the individual observations, and 
the mean obtained by a linear fit to all the observations 
(Sect. The data was scaled in time according to the 
solutions from Sect.0 The errors were calculated for each 
point directly from all contributing data, weighted accord- 
ing to the number of points (after time scaling) that con- 
tribute. 

Grand structure functions were calculated for all the 
data, and over subsets of interest, including in and across 
the direction of the elongation found above, and the indi- 
vidual years. We did not find any very significant differ- 
ence in the structure functions of these sub-sets. We show 
in Fig. ll3l the grand structure function obtained using the 
entire data set. 

The slope of the structure function before reaching sat- 
uration is 1.82±0.05 (i.e.{D'^{T)) cx t^-^'^). This should be 
compared with {{D^{t)) ~cx r) for classical IDV sources 




Fig. 11. The velocity vectors on the first of the month 
for the best fit model, using the constraints from Paper 
II. The velocity vectors are the velocity of the Earth 
through the scintillation pattern (f (Earth) — ^^(piasma)), as 
seen from the source. The dashed-dotted line shows the 
projected plasma velocity (?^(piasma) ) ■ In order to illustrate 
the combined effects of the velocity and the anisotropy we 
also present, on the same plot, the position angle of the 
anisotropy as seen from the source (solid line). 
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Fig. 12. The data adjusted by the timescales calculated from the best fit model. The timescale of the fluctuations is 
now constant. The apparent change in depth of the modulations is an effect of an increase in source brightness: the 
rms of the intensity fluctuations about the mean does not change significantly over the years (see Fig. inj. 



found by ISimonett i et all lIlQSfil) . iBlandford et alJ l)l986j) 

interpreted their result as due to scattering in a thick 
screen, which they showed was able to reproduce the 

cx T result, unlike thin scattering screens which give 

cx t2. 

On the basis of the slope of the structure function, 
it is therefore likely that the scattering in J1819+3845 
is predominantly in a relatively thin screen (thickness ^ 
distance). An inaccurate model for the annual timescale 
variations would decrease the slope of (D), but in the in- 
dividual runs we did not observe structure function slopes 
greater that 1.9. We also note that the grand structure 
function overshoots its saturation value. We consider the 
origin of these features. 



The D 



,1.82 



relation could potentially be explained 



if the screen, whilst thin, was allowed some finite thick- 
ness. We note that D (x for /3 = 4, (where (3 is the 
index of the power spectrum of the density fluctuations) 
and shallower power spectra (e.g. Kolmogorov) give rise 
to structure functions which rise more rapidly. Therefore 



shallo wer power spectra c annot account for the observa- 
tions feoma,Tii et ablllflSd Fig. 7). 



An anisotropic scatterer could produce a s tructure 
function which overshoots its saturation value jRickettl 
■200 LI . In this case we should expect this to clearly show 
up as a difference in the structure functions when mov- 
ing across, and along the direction of anisotropy. We see 
the overshoot when considering the entire dataset, and 
no strong evidence for a difference between the structure 
functions when sampling the scintles perpendicular and 
parallel to the elongation. We conclude that preliminary 
analysis does not indicate plasma anisotropy, however this 
result is tentative because of the small number of scin- 
tles which contribute to the GSF parallel to the elonga- 
tion. Source elongation/structure and scattering plasma 
anisotropy may be distinguishable by more sophisticated 
means (Rickett et al, in prep). 
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Fig. 13. Grand structure function, containing all the data, 
resampled according to the best fit from Sect. El The 
dashed line corresponds to the saturation value calculated 
from the modulation index calculated from all (resampled) 
data. 



8. The stability and structure of the source 

8. 1 . Source flux density changes 

The source has apparently undergone a change in mean 
flux density during the two years of monitoring. Fig. 
shows the mean flux densities, with the error on them 
calculated as in Sect. 4.1. Fitting a linear rate of change 
(using a weighted least squares flt) gives a rate of change 
of intensity as ~25% per year. A line with zero intensity 
change is a much poorer fit to the data. The change in 
flux density is also clearly seen in the distribution of ob- 
served flux densities (Fig. 114b). The mean flux densities 
per year (data now sampled according to the timescale 
flt) are 137 mJy (1999.8), 172mJy (2000.5), and 215mJy 
(2001.3). 

It is possible that the source became brighter at 
5 GHz in outbursts at certain epochs: in particular around 
February 2001 (MJD 51900). However, because of the 
extremity of the variations and the number of scintles ob- 
served per run, we cannot distinguish this possibility from 
a steady brightening. 

8.2. Source structure, and (lack of) changes 

We have no evidence for structure on arcsecond or ar- 
cminute scales in this source. From the radio spectrum 
(Paper I) we do not expect any such structures to con- 
tribute more than a few mJy. We have looked for evidence 
of extended structure by mapping several 12 hour runs. 
We removed J1819-f 3845 by amplitude self-calibrating the 
whole run on the mean flux density, and then looked for 
sources left in the field (which would now themselves ap- 
pear variable). There are a few faint compact sources 
more than 1 arcminute away, which have been seen also 
at 1.4 GHz. We find no evidence for any extended struc- 
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Fig. 14. (a)The mean flux density of source during the 
monitoring campaign. The solid line represents the best 
linear fit to the data; the dotted line the weighted 
mean.(b)The distribution of the observed flux densities in 
1999, 2000 and 2001. To avoid over-representation of cer- 
tain values, data from the slow periods has been sampled 
according to the timescales calculated previously. 

ture associated with J1819-f3845 to the level of 0.1 mJy 
for structure with scales 5-10 arcsec. 

From Sect. we concluded that the scintle had an 
elongated structure, but that this could be due to screen 
anisotropy and/or source structure. If the anisotropy is 
due to source structure it has been very stable, and by 
comparison of results from two year periods, we conclude 
that it has changed position angle by less than 10° in two 
years. 

Further information on the source structure can be 
gained from the probability distribution function (PDF) 
of its flux density. A point-like source has a very skew 
probability distribution function of the intensities, and a 
very jagged appearance of the time series. This is both ex- 
pected in theory and seen in observations of pulsars (e.g. 
iRickett et al.i.200()|) . The smoothness of the time series 
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- the quasi-oscillatory appearance - indicates that the 
source is resolved on the scale of a Fresnel zone at the 
scattering screen. 

There is also an indication of asymmetry in the source 
brightness. The falling edges of the flux density varia- 
tions exhibit a more rapid change than the rising edges 
l)de Bruvn fc Dennett-ThorDell200lt) . Initially a puzzling 
result, because there was no evidence of a change in sense 
over the year, this result can be explained with the scat- 
tering plasma velocity. The scintillation pattern is always 
cut in the same sense, and the Western end of the source 
is brighter than the Eastern. In a future paper we will 
return in more detail to the scintle asymmetry. 

We also use the PDF to look for changes in source 
structure over the observation period. We present the PDF 
for the entire data fFig. 114b'). where the data has been 
sampled according to the theoretical timescale from the 
model fit in The increase in brightness over the years 
is obvious. Normalizing on the mean for each year (or the 
linear change calculated above) , and further increasing the 
sampling rate to produce independent points (i.e. only one 
point in the calculated timescale is sampled) allows us to 
calculate statistical measures of difference. Years 1999 and 
2001 are statistically indistinguishable by the KS test (> r 
99% certainty). However year 2000 has a 20% probability 
(KS test) that the distribution is similar to the other two 
purely by chance. This is primarily due to the higher dis- 
persion about the mean (cr = 0.38, 0.44, 0.39 in 1999, 2000, 
& 2001 respectively). From this preliminary investigation 
we see no evidence for change in structure. 

In Paper I we considered the source to be composed 
of a scintillating component and a non-scintillating base 
component. However, the source may be a few times larger 
than the Fresnel scale, and all scintillating. This would 
also result in modulations in flux density that do not drop 
to zero, giving an appearance of a 'base level'. As previ- 
ously noted, the source is expected to be resolved in order 
to explain the smoothness of the variations. The 50% in- 
crease in source brightness with no concomitant change in 
modulation index, and indeed the lack of difference in the 
normalized flux density PDFs over the years, is difficult to 
explain in the 'baselevel plus small component' model. It 
would require that both the non-scintillating and the scin- 
tillating components increased in flux density by a similar 
amount, which is somewhat contrived, whereas this sit- 
uation is a natural result of a resolved source that gets 
brighter. For this reason we now favour the picture of a 
resolved source, a few times the Fresnel size. 

We then can interpret the modulation index in order 
to obtain the source size relative to the scattering zone, 
which we will return to in Sect. El 

8.3. Expansion rate 

From the modulation index alone (Fig.lH)), which shows no 
signs of changes over the two years, we can put an upper 
limit on the source expansion. As we have reason to believe 



that the source is not much smaller than the scattering 
disk (above), we can directly translate the lack of decrease 
in modulation index to a lack of detectable source expan- 
sion. Because the scintillation pattern is a convolution of 
the scintillation pattern produced by a point source and 
the source brightness distribution, an increase in source 
size would always act to decrease the observed modula- 
tions. In the strong sc attering regime , the modulation in- 
dex oc 1/ (source size) ()NaravaiJll992[) . We put a limit on 
the decrease in modulation index at 10%, and thereby cal- 
culate that the source has expanded by less than 10% over 
two years. 

9. Properties of the scatterer and source size 

In Paper I we used a simple formula based on the observed 
timescale, the frequency at which the modulations were 
maximum (~5 GHz), and an assumed velocity through the 
scintillation pattern to deduce the location and strength 
of the scattering screen and the size of the source, (note 
a typographic error in Paper I: calculated in that pa- 
per should read 0.017, not 0.17m~^°/^). We will return to 
this calculation again, now having ascertained the velocity 
through the scintillation pattern as a function of epoch. 

As the velocity of the plasma screen is of the order a 
few tens of km/s, and because of a lucky accident of the 
date of our first observations, the velocity used in Paper I 
was a reasonable figure, and therefore the measured veloc- 
ity of the plasma screen does not itself change significantly 
the apparent size of the source. However, we also need 
to take in account the effects of the source size and the 
anisotropy in the scintillation pattern. These introduce a 
number of significant changes. 

We reformulate the scattering results for a point source 
from Paper I, in terms of the desired quantities. These for- 
mulae are derived from Blandford & Narayan (1985) and 
Romani et al. (1986), using the point of maximum refrac- 
tive scattering (m=l), and assuming that the timescale is 
just the time taken to move through a scintle. 

0_4 - 1-^ t(J Km ^knis-^ '^hrs 
Lkpc = 1-9 10 tLs Km 

Here C^^IO'* C_4, 9 is the radius of the scattering disk, 
or the Fresnel disk at the critical frequency {v* — c/X*), 
and the distance to the equivalent screen is L, and we 
move with velocity v through the scintillation pattern with 
timescale t. These equations can only be applied directly 
to observations at the critical frequency for a point source, 
or to a source with a non-scintillating base component and 
a point- like scintillating component. 

As argued in Sect . 18 . 2l however . we consider the source 
to be resolved. The introduction of a source size means 
that the observables (m, t, v°) no longer reflect those of 
the refractive scintillations from a point source. We can- 
not therefore just drop the brightness temperature from 
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Fig. 15. The change of modulation index with frequency 
for source sizes larger than the scattering disk at the criti- 
cal frequency, v* . The solid line represents a point source, 
the dotted lines extended sources, with frequency depen- 
dencies as indicated. All extended sources are 3 times the 
Fresnel zone at the critical frequency {d*p). The points 
correspond to multifrequency observations, placed on the 
graph such that the 5GHz point is at v° for the case 

that calculated in Paper I by a factor 4 to compensate 
for a source size increase of a factor 2. In what follows 
we will turn the observed quantities into the equivalent 
point source values for use in the above equation, using 
the suffices * and ° to refer to quantities measured at the 
frequencies of maximum modulations for a point source, 
and for an extended source respectively. We refer to the 
frequency of maximum modulations for a point source, v* ^ 
as the 'critical frequency'. 

9.1. Effects of source size on m and t at maximum 
modulation 

A physical source would be expected to change its ex- 
tent with observing frequency. This has an effect on the 
characteristics of the depth of modulation as a function 
of frequency if the source is resolved {9s > 6'scatt). An ex- 
tended source does not exhibit maximum intensity vari- 
ations at the same frequency as a point source, but at a 
frequency which is also dependent on its change of size 
with frequency. This is because the observed frequency of 
maximum modulations for a point source occurs when the 
scattering disk is the same size as the Fresnel zone, but for 
an extended source this occurs when the scattering disk is 
equal to the source size. This is illustrated in Fig. [T51 

We are observing close to the frequency of maximum 
modulations. Simultaneous multifrequency observations 
from 10/2000 through 07/2001 confirm the effect reported 
in Paper I. The observations will be reported in a later 
paper, but the average modulation indices are plotted on 
Fig. El as it is critical to our analysis in this paper that 



the 5 GHz observations are close to the modulation max- 
imum. We can calculate the frequency of maximum mod- 
ulations that would be observed for a point source, if we 
assume that the source is larger than the scattering disk 
in the observed frequency regime, and that there is no 

non-scintillating component. 

We use the asymptotic theory l)Naravanll992l) . for both 
the strong and the weak regimes. Calculating the pre- 
dicted modulation indices for a resolved source, we find 
that for source size decreasing with frequency less rapidly 
than z^^^^/^, the observed frequency of maximum modu- 
lations is 

For illustration, the effect on the modulation index as a 
function of frequency for a number of source size frequency 
dependencies is illustrated in Fig. 1151 The frequency de- 
pendence of the source size does not, however, affect the 
above relation. 

Note that a non-scintillating component of the source 
would reduce the observed modulations without changing 
the observed modulation maximum from the critical fre- 
quency, but a resolved source results in a true critical fre- 
quency which is above the observed frequency of maximum 
modulation. This is because in this case, the maximum 
modulations do not occur at the transition from strong to 
weak scattering, but when the source size equals that of 
the scattering disk. Based on this simple calculation which 
assume the asymptotic theories are reasonable approxima- 
tions to reality, and assuming there is no non-scintillating 
component, and using the observed maximum in the in- 
tensity modulations (1^°) is around 5 GHz, we find that 
the critical frequency {v*) is around 25 GHz. 

We next consider the effect of source resolution on the 
observed timescale. Using the relation between the ob- 
served frequency of maximum modulation and the crit- 
ical frequency derived above, we find that the timescale 
observed at the maximum modulations is related to the 
timescale for a point source at the critical frequency, t*, 

by 

t° = (m°)-66/"i* 

Under the same assumptions as above, and using the ob- 
served m°=0.4 we find t* = t°/35. 

In order to relate the source size at the critical fre- 
quency to that at the frequency of observed modulation 
maximum, we must assume a dependency of source size on 
frequency. From multifreqency monitoring over the same 
period (forthcoming paper) we find that the spectral in- 
dices of the source around 5 GHz are Ofj'g ~ +1.0 and 
afg « +0.6. We therefore expect that the source will 
display a frequency dependence between that of an opti- 
cally thick source {a = +2.5, 9 constant) and a standard 
equipartition flat-spectrum source (a = 0, 9 oc We 
therefore consider a source which decreases in size as v~'^'^ 
(See also e.g. de Bruyn (1976)). 

In this case, and again for no non-scintillating compo- 
nent in the source, we find that the maximum modulation 
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index, m°, is related to the source size and the Fresnel 
scale at the critical frequency as 

Using the observed m°=0.4, we find 6*g = 1^0*p. Thus the 
source could be up to 15 times the size of the Fresnel scale, 
at the critical frequency, i.e. at 25 GHz. 

9.2. Effects of anisotropy 

Finally we must account for the anisotropy. Without con- 
sideration of the anisotropy we cannot turn the observed 
values into source and screen parameters, because of the 
order of magnitude change in observed timescale during 
a year. The interpretation of the anisotropy is the most 
difficult, and potentially the biggest source of error as the 
theory used is an isotropic scattering theory. In what fol- 
lows we simply take geometric means of relevant quanti- 
ties, and assume these would be the values appropriate 
to an isotropic theory. We consider two cases: that the 
anisotropy is due entirely to the medium, or entirely due 
to the source. It is possible that both the effects play a role: 
in this case the position angle of the anisotropy fitted to 
the data represents a position angle of the convolution of 
the two effects. 

If the anisotropy is due to the scattering medium, it is 
not clear how we should take this into account, as a full 
treatment of anisotropic scattering should be used, not 
the isotropic case on which we base our arguments. It is 
complicated by the fact that, as argued above, the source 
is likely to have a size comparable to the scattering scales. 
To make a step towards accounting for the anisotropy, we 
consider the situation when the Earth cuts through the 
anisotropy such that the scale is the geometric mean of 
the two scales (day 220, for best fit e= 14), when our 
velocity through the scintillation pattern is calculated to 
be 54km/s. 

9.3. Derived source and screen parameters 

First we consider that the medium is isotropic and the 
source is a resolved ellipsoid with no non-scintillating base 
component. Because of the conversion estimates we are 
making to treat the isotropic theory in a way to account 
for anisotropy, the results will also be those for the case in 
which the observed anisotropy is due to the medium, but 
the source itself is a partially resolved isotropic source. 

From the above considerations of the effect of source 
size we find t* = t°/35 — 0.057 hrs. We note that here 
we associate twice the timescales measured and presented 
in Table 01 with the time required in the formula, and 
the angular scales as corresponding to the diameter of the 
source. From the preceeding discussion, we associate t with 
t* and 9 with 9p in the formulae given at the beginning 
of Sect. 9. Using w=54km/s, m°=0.4, iy°=5GHz in these 
formulae, we find 9p= 25/ias. At 5 GHz the source is then 
0^ X 15 X = 860^as. The source has a mean 



fiux density of 200 mJy, resulting in a brightness temper- 
ature Tb = 2S\^/'Kk9l = 1.3 X 10^° K. In this case the 
screen is extremely strong (C^ ~ 165m~^°/^), and nearby 
('^ 1 pc) , which results in the larger angle subtended by 
the scattering disks. This in turn considerably lowers the 
brightness temperature from that calculated in Paper I. 

We now consider the case of the maximum possible 
flux density in a non-scintillating component. Although 
the modulation index is 0.4, we know that >40% of the 
source is scintillating, by the regular minima at ~ 60 mJy. 
We therefore consider, in a similar manner to Paper I a 
non-scintillating base-level determined by this value. From 
Fig. ll4r a) we take 140 mJy in the scintillating component. 
However our approach here differs from Paper I, not just 
due to our revised velocity, and calculated anisotropy, but 
because we take into account that this scintillating com- 
ponent must also be somewhat resolved. The reason for 
this is that the modulation index is not 0.7, as would 
be expected from scintillation of a 140 mJy point source 
at the critical frequency, on a non-scintillating compo- 
nent of 60 mJy. The observed modulation index of the 
smaller component, i.e. if we remove the base compo- 
nent, is m°' — m°/a;, where x is the fraction of the total 
flux density in the scintillating component. We use m°' in 
place of m° in the formulae given in the preceding sec- 
tions. Again using v=54km/s, t°=2hrs, and j/°=5GHz, 
we obtain at 5 GHz a source size of 100/ias and brightness 
temperature of 7 x 10^"'^ K. The scattering screen is 12 pc 
away with = 0.5m"20/3^ 

We have presented the boundary scenarios: the reality 
may lie somewhere in between. Fig. 1161 shows the two ex- 
treme cases discussed above {x=l and 0.7), as well as the 
intermediary cases. We note that the size quoted is the 
effective diameter, so for an axial ratio of 6:1, the source 
would be a/6 « 2.5 times this in the direction of the elon- 
gation, and ~ 0.4 times this orthogonally to it. 

Our observed v° at around 5 GHz falls near to the 
critical frequency predicted from contribu tions to scatter- 
ing from the entire Galaxv IWalkeJ l)l998(l . but this must 
be coincidental because, aside from the reasons given in 
this section, the peculiar velocity rules out a contribu- 
tions from a very extended material, as does the slope of 
the structure function. We calculate a high critical fre- 
quency of between 13 and 25 GHz (0.7 < a; < 1). In this 
particular direction we see that any source brighter than 
a hundred mJy at 5 GHz would have to have a Tb far in 
excess of the Compton limit if it were to demonstrate an 
observed modulation maximum around 25 GHz. The co- 
incidence is probably also partly a reflection of the more 
frequent observing at GHz frequencies, where this source 
is most variable. 

We consider also the possibility of the source as mod- 
elled by two well-separated sources, each smaller than the 
scattering disk but separated by several times this. Such 
a situation is similar to some structures observed in VLBI 
images of core-jets and compact radio galaxies. Motion 
parallel to the axis of separation of the components would 
produce a conspicuous 'dip' in the rising portion of the 
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Fig. 16. The derived range of source and screen properties 
(upper and lower panel respectively), taking into account 
the effects of source resolution on the observed timescale 
and modulation index at the frequency of maximum varia- 
tion. The solid lines refer to the scale on the left, the dotted 
lines to that on the right. We require that the combined 
effects of the non-scintillating component and the size of 
the scintillating component result in the observed maxi- 
mum modulation index, and assume 6 oc v~^'^ . The values 
&t X — 0.4 correspond to a point-like (^5 < ^scatt) com- 
ponent on a non-scintillating base. However, the left-hand 
side (shaded) is not compatible with other considerations 
of the data: the PDFs and the original lightcurves show 
this is not compatible with the deep minima: > 70% of 
the flux density must be in a scintillating component. 

structure function ijSpander et al.lll993(l . The nature of 
the observed structure functions for J1819-I-3845 indicates 
this is unlikely, as the fit indicate most observations as well 
as most scintles are observed during motion suitable to see 
the development of these two timescales. There is no such 
evidence for this: the 'overshoot' of the structure func- 
tion cannot be explained in these terms, as the structure 
function reaches saturation values before turning over. We 
therefore rule out a model of discrete and well-separated 
source components. 

Finally we note that, as well as not using anisotropic 
scattering theory, all these calculations are based on 
asymptotic theory which is not strictly applicable in this 
regime. Further work testing the theory with simulations 
near the critical frequency is required, as well as the ef- 
fects on the observables of partially and slightly resolved 
sources {63 ~ ^scatt)- Only more careful theoretical work 
will enable the scattering material and source size to be 
properly constrained. At present, even for a source with 
very well defined characteristics (a strong limit on 'non- 
scintillating' components, and a known velocity of the 
scattering screen) the constraints on source size are not 
tight. 

It can be seen that the brightness temperature can be 
below the Compton limit, but only if the scatterer is very 



strong and only a few parsecs from Earth. Regardless of 
the scattering theory, we note that if the source has a 
Ts < 5 X lO^^K, the source must be > 140^as in di- 
ameter. The detection of anisotropy, whether in scatterer 
or source, is a robust result which does not depend on 
theoretical details. The anisotropy required by the annual 
modulation fits has an axial ratio of > 6:1, and if this is due 
to source structure, the longest axis must be > 0.34 mas 
at 5 GHz. This is resolvable by space VLBI or ground- 
based VLBI at higher frequencies. Global VLBI observa- 
tions at 22 GHz (Dennett-Thorpe, Gurvits et al. in prep) 
show J1819-I-3845 to be <0.4mas at this frequency. If we 
assume that the source at 22 GHz closely related to that 
at 5 GHz, in particular the (not unreasonable) assump- 
tion that there are no milliarcsecond components which 
are present at 5 GHz and not at 22 GHz, we determine an 
upper limit on the source size at 5 GHz. Using the relation 
source size cx z^^*^'^, we find the major axis of the source 
at 5 GHz is < 8 mas, which is an equivalent diameter for 
axial ratio 6:1 of <0.32mas. From this we favour a source 
which exceeds 5 x lO^-'^K at 5 GHz. Combined with the pre- 
ceeding discussion this would put the screen at a distance 
between 4 and 12 pc. 

10. Concluding remarks 

We have presented observations over a two year period 
for the quasar J1819+3845. Notable results of this study 
were: 

(i) The source continues to vary throughout the two 
years. Although the timescale of these variation changes 
by a factor of ten or more, the strength of the modulation 
does not appear to change. 

(ii) The change in timescale has a periodicity of a year, 
and can be explained by a velocity of the plasma screen 
and an additional component, which is due to source struc- 
ture or anisotropic turbulence in the scattering plasma. 
These observations do not allow us to uniquely determine 
all four parameters, although Vdoc of the scattering plasma 
and the position angle of the anisotropy are well con- 
strained. Combining these observations with th e results 
of ou r previous experiment ( Dennett-Thorpe & d e Bruvnl 
I2OOII , Paper II) where we measured the velocity of the 
plasma screen directly, we constrain the axial ratio of the 
source, or the turbulent cells to be > 6:1 and oriented at 
PA=83°±4°. 

(iii) The source is resolved, as can be seen by the 
smoothness of the modulations. From the depth of the 
modulations, we deduce that the source at 5 GHz is at 
most a few times the scattering disk at the scattering 
screen. 

(iv) Based on the interpretation of the structure func- 
tion, we determine that the plasma responsible for scat- 
tering is a thin screen. 

(v) Any source expansion has been < 10% over the two 
years of monitoring. No significant change in structure has 
been detected. 
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(vi) The source has been gettmg brighter at 5 GHz. 
The increase in flux density has been ~ 25%/year. 
Coupled with the lack of source expansion, this 50% in- 
crease in flux density, must be largely due to an increase 
in the brightness temperature. 

Interstellar scintillation has already been proved to be 
the cause of the variations in J1819-f 3845 (Paper II), but 
this paper shows that monitoring observations over a year 
can also be a proof of scintillation. Since this project 
started this method has been applied to B 0917-1-624 
ifRickett et al.l200ltljauncev &: Macauartl200lj) . although 
using a very few points. Monitoring campaigns are 
more appropriate for the majority of IDV sources, as 
their timescales of modulation are not short enough to 
be able to conduct a two-telescope delay experiment. 
Furthermore, we demonstrate that the transverse veloc- 
ity of the scattering plasma can be calculated using this 
method: to our knowledge this is the first such measure- 
ment. 

Using the determined plasma velocity we obtain 
an equivalent angular source size of between 100 and 
900microarcseconds, and corresponding brightness tem- 
peratures in the range 10^° and lO^^K. In this calcula- 
tion we consider the effects of the spatial distribution of 
the scintillating flux - how much of the source could be 
smaller than the scattering disk. This consideration is the 
main reason for a decrease in T^ from lO^^K calculated in 
Paper I, although we still favour a source T^ > 5 x lO^^K. 
The scattering material is very unusual: it is very strong 
and between 1 and 12 parsecs away. 

Modelling of the source extent in relation to the the 
scattering disk is a critical factor in determining both 
the intrinsic source size and the parameters of the scat- 
tering medium, and we are undertaking work on this. 
Although we have taken account of source resolution, we 
have not used a full and proper approach to deal with the 
anisotropy, and this also needs to be addressed. 

This work emphasises that consideration of the effects 
of source extension could have a very big impact on the 
parameters derived. In the meantime we must be wary of 
using models of scintillating and non-scintillating compo- 
nents to describe scintillating extragalactic radio sources, 
as the source sizes and brightness temperatures obtained 
in such a way are very poor estimates. 

The detection of the anisotropy is a robust result which 
does not depend on the details of the theory, but whether 
the anisotropy is due to the source itself, or to the scat- 
tering material remains an open question. 
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